The influence of bulk evolution models on heavy-quark phenomenology 
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We study the impact of different Quark-Gluon Plasma expansion scenarios in heavy-ion collisions 
on spectra and elliptic flow of heavy quarks. For identical heavy-quark transport coefficients rela- 
tivistic Langevin simulations with different expansion scenarios can lead to appreciable variations 
in the calculated suppression and elliptic flow of the heavy-quark spectra, by up to a factor of two. 
A cross comparison with two sets of transport coefficients supports these findings, illustrating the 
importance of realistic expansion models for quantitative evaluations of heavy-quark observables in 
heavy-ion collisions. It also turns out that differences in freeze-out prescriptions and Langevin real- 
izations play a significant role in these variations. Light-quark observables are essential in reducing 
the uncertainties associated with the bulk-matter evolution, even though uncertainties due to the 
freeze-out prescription persist. 

PACS numbers: 14.65. Dw, 25.75.Ld, 24.10.Nz, 24.85. -|-p 



I. INTRODUCTION 

One of the striking diseoveries of the heavy-ion pro- 
gram at the Relativistic Heavy Ion Collider (RHIC) is 
that the medium created in 200 AGeV Au-Au collisions 
behaves like a nearly perfect liquid [iHll- To further in- 
vestigate this new state of matter, usually referred to as 
the Quark Gluon Plasma (QGP), penetrating and well 
calibrated probes are essential to quantitatively deduce 
the effect of the medium on those probes. One of these 
probes are heavy quarks. The experimental investigation 
of heavy-quark probes at RHIC has been ongoing for sev- 
eral years now [l^ . Two rather unexpected observa- 
tions have emerged. Heavy mesons, despite their large 
mass, exhibit (a) an elliptic flow comparable to that of 
light mesons, implying collective motion of heavy quarks 
in the expanding medium, and, (b) a suppression at high 
transverse momentum (pt) similar to light mesons, im- 
plying a substantial energy loss of fast heavy quarks while 
traversing the medium. 

Several theoretical approaches have been put forward 
to describe heavy-quark energy loss, or, more generally, 
heavy-quark diffusion, in the QGP [TTl42^ . Two basic 
ingredients are required to perform quantitative calcu- 
lations for the modifications of the initial spectra. On 
the one hand, one needs a good knowledge of the mi- 
croscopic interactions of heavy quarks in the plasma, as 
encoded in their transport properties (drag and diffusion 
coefficients). On the other hand, one needs a realistic 
description of the expanding medium through which the 
heavy quarks propagate. The experimental spectra of 
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heavy hadrons (and their decay electrons) only reflect a 
combination of both components. This may lead to am- 
biguities in disentangling the two ingredients. 

A closer inspection of the theoretical modeling reveals 
differences in both macroscopic expansion scenarios and 
microscopic transport coefficients. A first comparison of 
these issues has been carried out in Ref. [1^. It is the 
purpose of the present article to separate both ingre- 
dients and study the influence of the medium descrip- 
tion, the freeze-out prescription and the Langevin re- 
alization on the heavy-quark spectra. Here, we focus 
on two previously used medium descriptions, namely the 
hydrodynamical model by Kolb and Heinz [2Q], and the 
more schematic elliptical fireball model by van Hees et 
aL 

The objective of this work is not to reproduce exper- 
imental data for heavy-flavor observables. Rather, we 
employ a common model for the elementary interaction 
but change the medium through which the heavy quarks 
propagate, or alternatively use a common description of 
the medium and change the elementary interaction, to 
elucidate the influence of both ingredients on the final 
spectra. 

Our article is organized as follows: In Sec. |ll] we de- 
scribe the two medium descriptions used in our compar- 
ison and also briefly recall the main in gred ients to each 
of the two transport calculations [l^ Il6l. [23l . H^l- In 
Sec. mil we calculate the spectra for different combina- 
tions of transport coefficients and medium descriptions. 
In this way we can separate the influence of the expan- 
sion scenario from that of the elementary interaction of 
the heavy quarks with the plasma constituents. For dcf- 
initeness, this comparison is conducted with the numeri- 
cal (Langev in) implementation used in ReL [l^l, but we 
to the Langevin implementation as adopted 
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in Ref. [lal- In Sec. IIVI we elaborate on how light-quark 
observables, like the elliptic flow of pions, are constrained 
within the two expansion scenarios and can help to dis- 
tinguish the latter. We conclude in Sec. |Vl 



II. MEDIUM DESCRIPTIONS AND 
TRANSPORT COEFFICIENTS 

Let us start by recalling the basic features of the 
medium expansion in Ref. [l6| , where an elliptical fireball 
has been used to model the evolution of the medium cre- 
ated in & = 7 fm Au-|-Au colhsions at y/s — 200 GeV/c, 
referred to as the vHR (van Hees/Rapp) medium in the 
following. It employs a QGP equation of state with 
Nf = 2.5 effective flavors which yields an entropy density 
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with T being the temperature. The total entropy (S* ~ 
4600 in Ay = 1.8 units of rapidity) is assumed to be 
time independent. The fireball volume is parametrized 
as a function of time according to 



V{t) = Tra{t)b{t){zo + ct), 
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The tuning of the parameters, 

ao = 5.562 fm, A = 0.55 c/fm, &o = 4.450 fm, 
B = 1.3c/fm, Woc^^ 0.52c, Aij = 0.122c , 



(4) 



will be discussed below. Since the fireball is approx- 
imated as homogeneous, the thermodynamic variables 
only depend on time but not on position. Therefore, 
for a given volume evolution and total entropy, the tem- 
perature, pressure and energy density can be calculated 
as a function of time as 



P = —{l6 + 10.5Nf)T^-BciGP, 



e = Ts-p = 3- —(16 + 10.5Nf)T^ + Bqgp, 



(5) 
(6) 



where the bag constant, -Bqgp ~ 356 MeV/fm'^, ensures 
the continuity of pressure through the mixed phase (as 
in the EOS-Q of Refs. [H, [l^). The phase transition 



is modeled by a standard mixed-phase construction at 
constant temperature, Tc = 180 MeV, with critical QGP 
and HG energy densities of e^*^^ = 2.25 GeV/fm'^ and 
e^^ = 0.82 GeV/fm^ (the latter following fro m a ther- 
mal hadron- resonance-gas model). The expansion pa- 
rameters in Eq. ^ have been determined to mimic the 
time evolution of the hydrodynamical model of Ref. 
but with final quark-momentum spectra and elliptic flow 
(^2) adjusted to results of coalescence-model fits to em- 
pirical pion and kaon spectra . Note that this absorbs 
the effects of a subsequent hadronic evolution. The effec- 
tive quark mass has been set to niq — 0.3 GeV, and the 
freeze-out prescription has been chosen consistent with 
the equilibrium limit (i.e., long-time limit) of the post- 
point Ito realization of the Langevin simulation for the 
heavy quarks within the vHR model. 
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(cf. AppendixElfor details). This function is reminiscent 
of the one obtained within the so-called Milekhin freeze- 
out prescription [2^ [sO] and will therefore be referred to 
as "Milekhin-like" in the following. In Eq. ([7]), u{x) is 
the four- velocity vector field of the medium. In vHR the 
transverse fiow field, v±(t,x), is constructed with help 
of confocal elliptical coordinates in the transverse (x, y) 
plane, i.e.. 



r T 
v= I — Vh{t)cosv, — Va{t) sin v,0 



(8) 



where Va{t) = d{t), Vb{t) = b{t), and 

r = {x,y) = (0.6ao sinhw cosw, 0.6ao coshit sinw) 
■tb = {b{t) cos V, a{t) sin v) . (9) 

The expansion parameters in Eq. ([3]) produce final light- 
quark spectra with an average surface-flow velocity of 
Vbitmix) — 0.55 and total elliptic flow of V2 — 5.6%, at 
the end of the mixed phase. Note that these parameters 
are designed to reproduce hadron observables, not neces- 
sarily the result of a hydrodynamic model at the end of 
the mixed phase. 

The second model to describe the expansion of the 
medium is (2-|-l)D hydrodynamical calculation by Kolb 
and Heinz (labeled KH) [2al. The equations underlying 
ideal hydrodynamics are the conservation of energy and 
momentum and appropriate currents (e.g., baryon num- 
ber), 



9^T'^'^(x) = 0, d^rix)^0. 



(10) 



The equation of state is very similar to vHR, consisting 
of two parts. The hadronic phase is described as a gas 
of non-interacting hadronic resonances, containing essen- 
tially the same resonances as in the vHR flreball model 
but using a lower critical temperature, Tc = 164 MeV. 
Above the phase transition, in the QGP phase, the sys- 
tem is modeled as a non-interacting gas of u, d, s quarks 
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and gluons, with an external bag pressure (of a similar 
value as the one used in vHR) . The temperature and en- 
ergy density are local quantities in this approach. For 
more information on the modeling of the swtem with 
this approach the reader is referred to Ref. [23 . 

The parameters of both medium descriptions have 
been fixed to describe Au+Au collisions at ^/s = 
200 ^GeV at an impact parameter of 6 = 7 fm. Note 
that the hydrodynamical model describes the hadron- 
Pt spectra and V2 at thermal freeze-out in the hadronic 
phase, while the vHR fireball has been tuned to repro- 
duce quark-momentum distributions at the end of the 
mixed phase such as to reproduce hadronic spectra within 
a coalescence model for hadronization [28j . Thus, by con- 
struction, the elliptic flow, W2, in the vHR scenario is 
larger than in the KH scenario. 

In Fig. [T] we display the temperature-time profiles ex- 
tracted from KH hydrodynamics and the vHR fireball. 
In the hydrodynamical simulation the temperature is lo- 
cal, and for this comparison we take the temperature 
in the central cell. The temperature in the elliptic fire- 
ball depends only on its volume and is equal at all 
points in space. The vHR medium prescription starts 
at To = 0.33 fm/c and stops at roughly 5 fm/c with the 
end of the mixed phase. For the KH prescription the 
medium evolution lasts longer, starting at tq = 0.6 fm/c 
and ending the mixed phase at '--^ 7 fm/c, followed by a 
hadronic phase. The impact of the hadronic phase on the 
heavy-quark energy loss has been neglected when using 
the KH medium. 
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FIG. 1. (color online) Temperature profiles for KH hydrody- 
namics (central cell) and for the vHR elliptical fireball. One 
observes different lifetimes of the QGP and mixed phases for 
both approaches. Additionally, the hadronic phase of the KH 
description is depicted but heavy meson interactions are ne- 
glected. 

Next we turn to the different descriptions of heavy- 
quark diffusion. In Refs. [23l . [23 | a pQCD calculation with 
running coupling has been used, whereas in Ref. [l6j ma- 
trix elements based on an effective resonance model fl^ 
have been employed. For this study, the interactions are 
implemented through the use of a Langevin-transport 
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FIG. 2. (color online) Drag coefficients for pQCD with run- 
ning Qs (solid line) and for resonance-f pQCD interactions 
(dashed line), as a function of momentum for a temperature 
of 180 MeV. A clear difference between the two approaches is 
observed. 



model (cf. Appendix |^ for more details) characterized 
by ensembles of random momentum kicks. 



dp,- 



T{t,p)pjdt + VdiCjfc(i,p + ^dp)wk , (11) 



which shift the heavy-quark momentum at each time 
step. At] the first term, ^T{t,'p)pj , is a friction term, Wk 
are Gaussian-distributed random variables and the Cjk 
are related to the diffusion tensor, B (see appendix). The 
choice of ^ e {0, 1/2, 1} characterizes realizations of the 
Langevin process, known as prc-point Ito, Stratonovich 
and post-point Ito (or Hanggi-Klimontovich) relations, 
respectively. For a particle of mass, M, which is large 
compared to the temperature of the ambient medium, 
one has T ^ A + 0{T/M) « A independent of the re- 
alization, where A is the drag coefficient responsible for 
the energy loss of particle, defined from the microscopic 
interaction by 



^ 2pE J (27r)32fc J (27r)32fc' J (27r)32£;' "'^^^ 
X (27r)4<5W(p+fc-p'-fc')-J-^|A^.|' [p-p') . 



(12) 



Figure [2] depicts the drag coefficients computed in terms 
of the scattering-matrix elements of the underlying mi- 
croscopic model for the heavy-quark interaction as a 
function of three-momentum at a fixed temperature of 
T = 180 MeV for the resonance model of Ref. [l2| and 
for the pQGD model of Refs. (dubbed TAMU and 

SUBATECH coefficients, respectively). One finds sizable 
differences by more than a factor of two. 
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III. SYSTEMATIC COMPARISON OF 
HEAVY-QUARK SPECTRA 

In this section we analyze heavy-quark spectra and el- 
hptic flow by combining both medium descriptions with 
both transport coefficients, as alluded to in the previous 
section. In particular, we also study the time evolution 
of the elliptic flow. This requires the identification of a 
suitable variable characterizing the time evolution. In 
the vHR medium the energy density and the evolution 
time are uniquely correlated, but such a definitive rela- 
tion does not exist in the hydrodynamic approach. In 
the KH medium we terminate the interactions of a heavy 
quark in a fiuid cell as soon as its energy density falls 
below the freeze-out energy density, efo , and evaluate the 
heavy-quark spectra on the pertinent hyper-surface. 

In Fig. [3] we show the total V2 for c-quarks, defined by 



/ prdprdip cos(2(/3) 
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(13) 
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as a function of different freeze-out energy densities, efo, 
down to the last point, where still a fraction of QGP ex- 
ists, Cfo = 0.45GeV/fm3 in KH and efo = 0.82GeV/fm3 
in vHR. All curves are calculated with the same micro- 
scopic model for the heavy-quark transport coefficients 
(pQCD with running coupling) [2^, 12^1. but with differ- 
ent medium evolutions and Langevin realizations. The 
full red line (with diamonds) depicts the KH hydro evo- 
lution with the pre-point Langevin, the dashed line (with 
triangles) the equivalent result for the vHR fireball, and 
the full blue line (with circles) the vHR medium with a 
post-point Langevin; the latter corresponds to the results 
of Refs. [ii. 
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FIG. 3. (color online) Total V2 of charm quarks as a function 
of freeze-out energy density for two different medium descrip- 
tions with identical (pre-point) Langevin realization (triangles 
and diamonds) and for the vHR evolution with post-point 
Langevin (circles). 



We first note that the precise realization of the 
Langevin process has a moderate influence on the out- 
come of the calculation in the heavy-quark sector. For 
the two different medium descriptions, the discrepancies 
are larger, up to a factor of 2 at a given energy density. 
This suggests that both media carry different momen- 
tum eccentricities along their evolutior0, which, in turn, 
are transferred to the heavy-quark motion. We already 
pointed out that the vHR evolution has been tuned to 
the empirical V2 of light quarks at the end of the mixed 
phase. We will return to this question in Sec. IIVI below: 
apparently, the moderate difference between 5.6% (vHR) 
and 4.8% (KH hydro) is not the main cause for the effect 
seen in Fig. [3] 

In the following, we elaborate on how the differences 
in the medium expansion affect heavy-quark spectra by 
performing a comparison with different transport coeffi- 
cients and medium prescriptions. Specifically, we com- 
pare two basic features of heavy-quark spectra at RHIC, 
i.e., the nuclear suppression factor, Raa{pt), and the 
elliptic flow, U2(pt), by studying four setups: 

I. KH medium / SUBATECH coefficients 

II. KH medium / TAMU coefficients 

HI. vHR medium / SUBATECH coefficients 

IV. vHR medium / TAMU coefficients 

The scenario described as "SUBATECH coefficients" are 
those published in Refs. [23l. [23|. whereas the drag and 
diffusion coefficients used in Refs. [H, [3 labeled 
as "TAMU coefficients". The scenarios I and IV de- 
scribe the published data on the nuclear suppression fac- 
tor, Raa, of semileptonic electrons fairly well, despite of 
different assumptions on the medium and the transport 
coefficients. To disentangle the effects of the two compo- 
nents, we swap the medium description with the descrip- 
tion of the heavy-quark interaction (scenarios II and HI). 
In this way we can cross-compare the results and better 
identify changes due to the medium or the diffusion mech- 
anism. Since the precise realization of the Langevin pro- 
cess has rather little impact in the heavy-quark sector, we 
have chosen to proceed with the pre-point (post-point) 
prescription for the KH (vHR) medium; the motivation 
for this choice will be justified in Sec. IIVI 

Figure U] depicts the nuclear suppression factor of 
charm quarks at the end of the respective mixed phases 
as a function of their transverse momentum for the four 
scenarios mentioned above. Scenarios I and IV give 
rather similar results, within ca. 20%. When changing 
the medium description and leaving the diffusion mecha- 
nism the same, i.e., comparing scenarios I with HI, or II 



This feature, which we later on refer to as "intrinsic V2" of 
^ with the longitudinal coefficient imposed as Bj^ = TEA in the medium, could be further quantified by, e.g., evaluating the 

order to enforce the Einstein relation. anisotropy, ep, of the energy- momentum tensor. 
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FIG. 4. (color online) Nuclear modification factor, Raa, of 
charm quarks at RHIC energies for the four different scenar- 
ios. The lines are identified as scenario I (full/thick red), 
scenario 11 (dashed/thick red), scenario 111 (full/thin blue) 
and scenario IV (dashed/thin blue). 



with IV, one observes a difference of 30-50% at high pr 
in opposite directions, indicating that the vHR medium 
induces more "stopping" than the KH medium. As ex- 
pected, the smaller friction coefficients of Refs. [l^, fl6| 
cause a smaller energy loss than the pQCD-|-running-as 
model. The maximal deviation of roughly a factor of 
2 occurs when combining the large coefficients with the 
"more stopping" medium (scenario III) compared to the 
small coefficients plus "less stopping" medium (scenario 
11). 

Similar features are found when comparing the ellip- 
tic flow of charm quarks, V2{pt), see Fig. [5] One ob- 
serves again that scenarios I and IV compare reasonably 
well to each other, on the level of ca. 20-30%. However, 
when interchanging the medium description the elliptic 
flow either increases by more than a factor of 2 (scenario 
HI compared to I) or decreases by roughly a factor of 2 
(scenario II compared to IV) at pr — 2-4 GeV. The dif- 
ferences at intermediate pt are thus stronger for V2 than 
for Raa- 



IV. ELLIPTIC FLOW OF LIGHT CHARGED 
PARTICLES 

As mentioned above, the calculated Raa and V2 val- 
ues of charm quarks cause ambiguities in quantitatively 
disentangling the effects of the expansion scenario from 
that of the microscopic interaction of heavy quarks in 
the QGP. It is therefore important to return to light- 
quark observables to better distinguish between different 
scenarios. Light-hadron observables depend on the fl- 
nal state of the expansion scenario but arc independent 
of the interaction of heavy quarks with the plasma con- 
stituents. Thus they allow to scrutinize the description 
of the medium evolution if the freeze out description were 
unique, which is unfortunately not the case. A key ob- 
servable is the elliptic flow of pions or charged particles. 
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FIG. 5. (color online) Elliptic flow of charm quarks at RHIC 
energies. The red (blue) lines are computed with the KH 
hydro (vHR fireball) evolution, and the solid (dashed) corre- 
spond to using pQCD+running-Qfa (resonance) model trans- 
port coefficients. 



In the following we study in more detail how our two 
bulk evolution iTiodels have been adjusted to experimen- 
tal data. 

We adopt the standard definition for the differential 
elliptic flow of a particle with mass, m. 



V2{pt) 



(14) 



in terms of the single-particle momentum distribution 
functior0 



d^iV 



PTdpTdip 



Ed^N 
d^p 
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If the expanding medium is in local thermal equilibrium 
during the expansion the key question is, how to convert 
the fluid cells, characterized by a temperature and flow 
field, into a particle distribution. 

In the KH hydro calculations (as in iTiost other hy- 
drodynamical models) the Coopcr-Fryc (CF) prescrip- 
tion [3l| is employed to evaluate the momentum distri- 
butions of particles after freeze-out. It converts a thermal 
medium instantaneously into a momentum distribution 
given by 



Ed^N 
d^p 



da^ pf^f{p,T,u). 



with the Boltzmann - Jiittncr distribution, 



f{p,T,u) 



1 



(27r)= 



■ exp 



p ■ u 



(16) 



(17) 



^ The integration over rapidity, y, is performed to collect the flow 
of all light partons that are able to interact with the heavy quarks 
in the rapidity interval Aj/ = 1.8 inherent to the fireball. 
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where T and u are the temperature and four-velocity at 
thermal freeze-out below which no further interaction oc- 
curs, and (T^ is the hypcrsurface at constant T (or energy 
density). In the following, we will consider this prescrip- 
tion to analyze quark spectra during the evolution of the 
QGP as well. 

In the vHR fireball model, the Milekhin-like freeze- 
out prescription has been adopted. Since the medium is 
approximated as isotropic, the hypersurface corresponds 
to the entire ellipsoid volume defined by the condition 



y 



< 1 



(18) 



The four- velocity at time, t, and positi on, x±, is evalu- 
ated from Eq. According to Refs. [U [s^l one thus 

has 



Ed^N 



dyprdprdip d^p 



dV 



i2n) 



3 E /ml (p, T, u) , 



(19) 

with /ml defined in Eq. ([T])- In the following, we will 
also evaluate the fircball-U2 with the CP prescription, by 
taking /ml ^ / in Eq. 

In Fig. [6] we display the differential V2 for constituent 
quarks with a mass m = 300 MeV, emanating from the 
KH and vHR medium, using both the Cooper-Frye and 
the Milekhin-like descriptions for the latter. Since the 
heavy-quark evolutions are terminated in both media at 
the end of the mixed phase, we choose the correspond- 
ing hypersurfaces for comparison. It turns out that the 
KH-I-Coopcr-Fryc and vHr-f-Milckhin-like prescriptions 
are rather close up to pr — 700 MeV. This encompasses 
most of the bulk-particles in the medium (and most of 
the interactions of the heavy quarks occur with those soft 
partons). However, if one applies the CF freeze-out to the 
vHR fireball one finds that the parton V2 is systematically 
above the KH+CF and the vHR-l-ML medium, even at 
low pt (e.g., by almost a factor of 2 at pr — 700 MeV, 
where KH-I-CF and vHR-l-ML cross). This reiterates the 
evidence found in the context of the charm-quark V2 in 
Sec, mil that the 'i;2-content of the vHR medium is signif- 
icantly larger than that of the KH medium. 

It is furthermore instructive to examine the total par- 
ticle V2 according to Eq. (|13|) - averaged over rapidity as 
in Eq. (|15p -, which is particularly suitable to illuminate 
time (or energy-density) dependencies. In Fig. [7] we dis- 
play this quantity as a function of the freeze-out energy 
density for m = 300 MeV "partons" above the critical 
value (0.45 fm"^ and 0.82 fm"^ for KH and vHR, re- 
spectively) and for m = 140 MeV "pious" below. 

In the KH medium, for e < 1.4 GeV/fm^, the 
direct evaluation of V2 from the physical fields (e.g. 
2kh('''b, 2?^)) is hindered by numerical fluctuations inher- 
ent to Vrfo(a;^), where Tfo(a;^) is the freeze-out Bjorkcn 
time for a given position in transverse space. We have 
therefore resorted to a direct Monte-Carlo sampling of 
the freeze-out hypersurface, S, in order to evaluate this 
quantity in a more robust way. 
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FIG. 6. (color online) Elliptic flow of particles with mass, 
m = 300 MeV, as a function of transverse momentum at 
the end of the respective mixed phase for the KH medium 
with Cooper-Frye freeze-out (solid red line) and the vHR 
medium-|- Milekhin-like freeze-out (solid blue line). The 
dashed line is the result of applying Cooper-Frye freeze-out 
to the vHR medium. 



For the vHR fireball, there is no problem in perform- 
ing the calculation for W2(e) down to the end of the mixed 
phase (cfo = 0.82 GeV/fm'^). By construction in the orig- 
inal work [l^ the Milekhin-like freeze-out prescription 
results in an total V2 of close to 6%, adjusted to experi- 
ment in semicentral Au-Au collisions at RHIC. The same 
is true for the KH-hydro-f-Cooper-Frye freeze-out, where, 
however, only 3/4 of the value is reached at the end of the 
mixed phase while the remaining 1 /4 of the experimen- 
tal value develops in the hadronic phase. Applying the 
CF freeze-out to the fireball leads to significantly larger 
values of ca. 9% at the end of the mixed phase, which is 
« 60% (90%) larger than the vHR-t-ML (KH+CF) pre- 
scriptions. 

Pursuing the bulk evolution beyond the end of the 
mixed phase, one can compare the elliptic flow of directly 
produced pions, v^'^{tt^). In the case of the KH hydro, 
we have evaluated the 172 directly from the tt"*" spectra re- 
sulting from the hydrodynamical evolution. For the vHR 
medium, we have proceeded as described above for par- 
tons but simply taken m = 140 MeV and a Boltzmann 
distribution to resemble pions after the mixed phase. The 
resulting pion V2 at the end of the mixed phase turns out 
to be 4.5% (8%) for Milekhin-like (CF) freeze-out, which 
is smaller (larger) than for the KH medium by roughly 
25% (30%). 

The analysis of the bulk V2 summarized in Figs.|n]and[7] 
suggests an explanation for the origin of the discrepancies 
found in the charm-quark spectra analyzed in Sect. IIIII 
On the one hand, different frcezc-out (and hadronization) 
prescriptions have been applied in the light-quark sector 
for the two medium descriptions, both of which lead to 
good agreement with the empirical pion V2 ■ On the other 
hand, when the two medium evolutions are analyzed with 
the same freeze-out prescription, an appreciable discrep- 
ancy in the bulk V2 emerges. This "intrinsic" bulk V2 
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FIG. 7. (color online) Total V2 as a function of freeze-out 
energy density for constituent quarks (m — 300 MeV, open 
symbols) and pions (m — 140 MeV, full symbols). The light 
quarks are evaluated with CF freeze-out for the KH (dia- 
monds) and vHR (triangles) medium, as well as with the 
original Milekhin freeze-out for the vHR medium (circles). 
The colored boxes at the bottom sketch the three stages of 
the medium evolution (QGP, mixed and hadronic phase). The 
pions have been evaluated in the hadronic phase (using Boltz- 
mann statistics, no resonance feeddown) either at the end of 
the mixed phase (vHR medium) or throughout the hadronic 
evolution for the KH medium. 



appears to be the key quantity in communicating the 
momentum anisotropy to the heavy quarks propagating 
through the medium. Maybe somewhat surprisingly, the 
results in the heavy-quark sector, for "realistic" trans- 
port coefficients, are not very sensitive to the Langevin 
implementation (i.e., pre-point vs. post-point prescrip- 
tion), which, in turn, dictated the chosen freeze-out pre- 
scriptions in both media in the first place. Clearly, the 
theoretical issue remains to better understand the dis- 
crepancies in the different freeze-out descriptions in con- 
nection with the underlying Langevin implementation of 
heavy-quark diffusion, which of course should be consis- 
tent. 



V. CONCLUSIONS 

We have discussed the impact of the bulk-medium evo- 
lution on heavy-quark phenomenology in heavy-ion col- 
lisions by studying two expansion scenarios, which have 
been applied earlier to calculate the elliptic flow and nu- 
clear suppression factor of heavy-quark observables at 
RHIC. In both approaches the final results depend on 
the microscopic interaction of the heavy quarks with the 
QGP (as encoded in their transport coefficients) and on 
the expansion scenario which provides different intrinsic 
V2 values of the bulk medium. We separated these ingre- 
dients in order to better understand the influence of the 
medium descriptions on final spectra. By switching the 



medium and freeze-out descriptions one observes differ- 
ences of around 50% in the Raa and V2 of charm-quark 
spectra. This effect has been cross-checked with different 
models for the microscopic input for drag and diffusion 
coefficients. 

In the past, many efforts have concentrated on a better 
understanding of the microscopic interactions between 
heavy quarks and plasma constituents. Our study sug- 
gests that the influence of different expansion scenarios 
on the heavy-quark observables is comparable to that of 
different descriptions of the heavy-quark transport coef- 
ficients. In principle, light-quark observables can help 
to determine the expansion scenario but differences in 
freeze-out prescriptions and Langevin implementations 
induce significant uncertainties at present. Different 
freeze-out descriptions can indeed result in similar values 
for the light-meson V2 in models where the "intrinsic" 
elliptic flow of the hot medium differs appreciably. 
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Appendix A: Langevin simulations and freeze-out 
scenarios 

One of the difficulties in the use of relativistic Langevin 
simulations for heavy-quark diffusion processes is the de- 
pendence of the resulting phase-space distribution func- 
tion on the realization of the stochastic integral. The 
Langevin process is defined by the time step, 



dpj = -T{t,p)pjdt + VdiCjk{t,p + ^dp)wk , 

for the heavy-quark position and momentum coordinates 
with respect to the rest frame of the heat bath. The Wk (t) 
denote stochastically independent normally distributed 
random variables ( "white noise" ) , 

{wj{t)wkit')) =S{t-t')S,k ■ (A2) 

As elaborated in Ref. [2^, the Langevin process defined 
by Eq. (jAl[) is equivalent to the Fokker-Planck equation 



for the heavy-quark phase-space distribution function /, 
where one can identify the usual drag force Aj = pjT — 



^Cik^^f- as well as the diffusion tensor B = C ■ /2. 
Here, ^ G [0, 1] defines the realization of the stochastic in- 
tegral, incorporating the effects of the force fluctuations 
around the average drag or friction force, governed by 
the drag coefficient, F. In Ref. [l^, a pre-point Ito pre- 
scription has been adopted = 0), in which case one has 
identically T = A. Furthermore, the diffusion coefficients 
Bl and Bt have been adjusted to enforce the Einstein 
relation by solving Eq. (18) of Ref. [s^ and imposing the 
additional constraint that 



Bt{p) 
Bl{p) 



obrutc 

JDrp 

^brutc 



ip) 



(A4) 



where i^^rute S^rute ^^ic coefficients evaluated 
directly with help of the differential cross-section for 
the microscopic q/g + Q q' /g' + Q' processes. This 
prescription has been adopted in order to preserve the 
anisotropy observed in high-energy collisions. Once these 
assumptions were adopted for the case of a fluid at rest, it 
has been checked numerically that the equilibrium limit 
of the distribution in a moving fluid is compatible with 



/BJlPiab) oc exp — j 



(A5) 



where u = u(t, x) is the four- velocity flow field of 
the background medium with respect to the laboratory 
frame. 

In Ref. [H, [ll] , following Ref. [11] , the Langevin real- 
ization has been chosen so that heavy quarks reach ther- 
mal equilibrium in the long-time limit with the temper- 
ature given by the surrounding medium, leading to the 
Boltzmann-Jiittner distribution, 



/oq(p) oc exp(- 



(A6) 



To avoid the evaluation of momentum derivatives of 
the diffusion coefficients, the post-point Ito realization 
has been adopted for the stochastic integral = 1). 
This allows to set F = ^ with the drag coefficient 



given by Eq. ([T^ . As for the pre-point Ito, the longi- 
tudinal drag coefficient has been enforced to obey the 
Einstein-dissipation-fluctuation relation, Bl = TEA, 
with E = \J + p^, which indeed leads to the equilib- 
rium limit, Eq. (|A6|) . independently of the speciflc mo- 
mentum dependence of the drag coefficient A j2H. For 
a flowing background medium, flrst the momentum co- 
ordinates have been Lorentz-boosted to the local heat- 
bath rest frame. After performing the time step (|Aip 
in this frame, the new momentum variables have been 
transformed back to laboratory-frame coordinates. This 
procedure leads to an equilibrium limit 



c ,^ N Plab • It "lab- 
/cq(Plab) OC e 



ab 



(A7) 



In terms of a Cooper-Frye freeze-out description, 
Eq. (jTB]) . this corresponds to the choice of a hyper- 
surface element da^ = d'^x w^, which is similar to 
the modified Milekhin-freeze-out description discussed in 
Refs. [l^H^I. In order to render the thermal-fireball de- 
scription of the bulk consistent with the freeze-out pre- 
scription implied by the equilibrium limit of the Langevin 
realization, the elliptic fireball described in Sec. |ll] has 
been adjusted using Eq. (jA7p as the local-equilibrium 
distribution of the light quarks. This leads to a light- 
quark bulk elliptic flow of V2 — 5.6%, as shown by the 
endpoint of the Milekhin-like curve in Fig. [T] This value 
is compatible with the experimental vi for light hadrons 
when applying the coalescence model for hadronization of 
Ref. (281 (the same model has been subsequently used in 
Refs. |l6l. [l9| to convert the heavy-quark spectra from 
the Langevin simulations into D- and i?-meson spec- 
tra). As pointed out in the text, if instead the standard 
CF constant-lab time freeze-out prescription is used, for 
which do" = d'^.T(l, 0, 0, 0), the same fireball evolution 
leads to a higher bulk elliptic flow of W2 — 9.2% for light 
quarks of mass = 300 MeV. This finding reiterates 
a main point of the present investigation: Conclusions 
about the microscopic dynamics (i.e, the transport coef- 
ficients of heavy quarks in the QGP) from experimental 
heavy-quark observables, like Raa and Vi of single elec- 
trons at RHIC, depend on the description of the expan- 
sion of the QGP. This expansion has to be understood 
from the light-hadron data. 
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